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ABSTRACT 

I find evidence for clustering in age of well-dated impact craters over the last 500 Myr. 
At least nine impact episodes are identified, with durations whose upper limits are set 
by the dating accuracy of the craters. Their amplitudes and frequency are inconsis¬ 
tent with an origin in asteroid breakups or Oort cloud disturbances, but are consistent 
with the arrival and disintegration in near-Earth orbits of rare, giant comets, mainly 
in transit from the Centaur population into the Jupiter family and Encke regions. 
About 1 in 10 Centaurs in Chiron-like orbits enter Earth-crossing epochs, usually 
repeatedly, each such epoch being generally of a few thousand years duration. On 
time-scales of geological interest, debris from their breakup may increase the mass 
of the near-Earth interplanetary environment by two or three orders of magnitude, 
yielding repeated episodes of bombardment and stratospheric dusting. I find a strong 
correlation between these bombardment episodes and major biostratigraphic and geo¬ 
logical boundaries, and propose that episodes of extinction are most effectively driven 
by prolonged encounters with meteoroid streams during bombardment episodes. Pos¬ 
sible mechanisms are discussed. 

Key words: comets: general - interplanetary medium - meteorites, meteors, mete¬ 
oroids - zodiacal dust 


1 INTRODUCTION 

The discovery of an excess of iridium at the Cretaceous- 
Tertiary boundary of 65 Myr BP, and its interpretation as 
the result of an asteroid impact (Smit & Hertogen 1980, Al¬ 
varez et al 1980), lent support to the hypothesis that large- 
body impacts are the prime cause of mass extinctions on 
Earth (Napier & Clube 1979, Raup 1991, Rampino 1994). 
However subsequent studies have not, with the still-debated 
exception of the Cretaceous-Tertiary event itself, confirmed 
this hypothesis. For example the Chesapeake (40 km) and 
Popigai (90 km) craters are amongst the largest known ter¬ 
restrial impact structures, dated to ~35.5 Myr, in the late 
Eocene, but are associated with only minor, regional extinc¬ 
tions of a few radiolarian species. The 80 km Manicougan 
crater, aged 214 Myr, is likewise associated with only mi¬ 
nor extinctions. Matsumoto & Kubotani (1996) found that 
the mass extinction and impact crater records are correlated 
with only 93-97% significance. 

Considerable attention has been paid to the question 
of periodicity in the impact cratering record, and whether, 
if such exists, it can be connected with Calactic cycles or 
mass extinction episodes; the issue remains controversial (cf 
Bailer-Jones 2009, 2011, Melott & Bambach 2013). Less at- 
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tention has been paid to the issue of episodes of bombard¬ 
ment, although several such have been noted and claims 
have been made of correlation with mass extinction (Keller 
et al 2012; McGhee 1996). Comet showers induced by stel¬ 
lar penetration of the inner Oort cloud, or asteroid showers 
deriving from the breakup of main belt asteroids, are inad¬ 
equate by about an order of magnitude to reproduce the 
sharpness, amplitude and frequency of the bombardment 
episodes (Napier & Asher 2009). Leakage from the main 
asteroid belt appears in any case to be inadequate by an or¬ 
der of magnitude to supply near-Earth object impactors say 
more than 3 or 4 km across (Menichella et al 1996, Minton 
& Malhotra 2010). It has been proposed that bombardment 
episodes are caused by the disintegration of giant comets in 
short-period. Earth-crossing orbits (Clube & Napier 1984). 
A more recent worldwide extinction of species at 12,800 BP, 
accompanied by a climate cooling which lasted for over 1000 
years, has likewise been attributed to an encounter with the 
debris from such a body (Napier 2010; Wittke et al 2013; 
Moore et al 2014), although this remains a controversial 
proposition (van Hoesel et al 2014). Models that neglect 
the injection of large comets into the near-Earth environ¬ 
ment (e.g Chapman 2004, Boslough et al 2013) fail to yield 
large-scale catastrophic disruption of this sort in the recent 
past, but it remains to be seen how often such comets ar¬ 
rive and what effect they may have on the biosphere over 
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geological as well as human timescales; this is the subject of 
the present paper. 


2 TRANSFER INTO SHORT-PERIOD ORBITS 

Recent years have seen the discovery of large populations 
of apparently icy bodies, some in unstable orbits, at the 
edge of the planetary system. Hahn & Bailey (1990) have 
demonstrated that 95P/Chiron, the first Centaur to be dis¬ 
covered, is in a chaotic orbit and may have entered a short- 
period, Earth-crossing configuration several times over the 
past ~10® yr. From Spitzer observations its diameter is esti¬ 
mated to be 233±14 km (Stansberry et al 2008). The poten¬ 
tial for icy bodies of 100 km diameter upwards to enter short 
period orbits with perihelia q <1 au has been demonstrated 
both theoretically (Hahn & Bailey loc. cit.) and observa- 
tionally. Sekanina and colleagues (2000, 2003, 2004, 2007) 
have shown that the Kreutz group of sungrazing comets is 
derived from a complex hierarchy of fragmentations, start¬ 
ing with the disintegration of a progenitor comet of probable 
diameter ~100 km about 1700 years ago. Likewise an exten¬ 
sive study of the Taurid complex supports the view that it 
originated from the breakup of an exceptionally large comet 
in the distant past (Whipple 1940, Sekanina 1972, Clube & 
Napier 1984, Asher et al 1994). The currently observed sys¬ 
tem may be at least 20,000 years old (Steel et al 1991), but 
an injection in a much more distant past is also possible. 

A 100 km comet striking the Earth would carry ~1000 
times the energy involved in the creation of the 150 km 
Chicxulub crater and would presumably remove the surface 
biosphere; impact crater statistics show this to have been an 
unlikely event over the 600 Myr timescale of the Phanero- 
zoic. However, disintegration seems to be a normal if not 
predominant mode of comet loss. Given that a Chiron-sized 
body has a mass lO^^-lO^^ gm, as against < 10^® gm for the 
entire current near-Earth asteroid system, there is the poten¬ 
tial for a substantially enhanced hazard during the breakup 
of such a body in the inner planetary system. Such a breakup 
might involve multiple bombardment episodes and substan¬ 
tial dust input (Clube & Napier 1984), leading perhaps to a 
qualitatively different astronomical scenario for mass extinc¬ 
tions from that of a single, 10 km asteroid impact (Alvarez 
et al 1980). 

2.1 Numerical integrations 

Di Sisto & Brunini (2007) modelled the Centaur popula¬ 
tion on the assumption that it derives from the scattered 
disc one, following the evolution of 1000 objects for 4.5 
Gyr and estimating the population by scaling from an as¬ 
sumed scattered disc population of 7.5 x 10® bodies with radii 
R >1 km (Fernandez et al 2004). They deduced a pop¬ 
ulation of 2.8x10® Centaurs, with an order of magnitude 
uncertainty. The mean dynamical lifetime Ld of Centaurs 
was found to decline strongly with perihelion distance, with 
Ld <10 Myr for objects with q <17 au, down to 7.6 Myr for 
bodies with 5.2< a <30 au and ~6.2Myr for Centaurs de¬ 
fined as Q <40 au and a <29 au. Horner et al (2004a, 2004b) 
obtained Ld ~2.7 Myr for a sample of Centaurs defined by 
the latter criteria. 

The hazard represented by the arrival of large comets is 
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Figure 1. Survival times of 86 out of 100 Chiron clones before 
ejection from the solar system. 
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Figure 2. Distribution of the time spent by 100 Chiron clones 
in orbits with q < 1 au. The 65 Earth-crossing epochs were due 
to 9 clones entering and re-entering these orbits. 


determined in part by their rate of arrival in the near-Earth 
environment rather than their mean lifetime as Centaurs, 
and this rate is dominated by that of the most dynamically 
unstable bodies at any given time. Horner et al (2004b), 
simulating the evolution of Chiron with 729 clones for up to 
3 Myr, found that over this time, more than half the clones 
became short-period comets at some point, and 85 became 
Earth-crossers. Similar behaviour occurred in other Centaur 
clones with similar orbits. 

I carried out orbital integrations of Chiron clones over 
1 Myr with the Mercury software developed by Chambers 
(1999), using the same input conditions as Horner et al 
(2004b). All the planets were incorporated. The hybrid part 
of the package was used (employing a symplectic integrator, 
but switching to Bulirsch-Stoer for close encounters). The 
accuracy parameter was set to 10“^® and initial timestep 
was 60 days; over a typical run this conserved energy to 
a fraction ~1.7 x 10~® and angular momentum to 10“^®. 
Of 100 clones, 86 were ejected within 1 Myr, their mean 
dynamical lifetimes before ejection being 0.35 Myr (Fig. [T]). 

© 2015 RAS, MNRAS 000, ?? - ?? 













Giant comets and mass extinctions 3 



0 100000 200000 


time (yr) 



time (yr) 

Figure 3. Top: Chiron clone to Earth-crosser. The upper curve 
represents the semimajor axis, the lower one is the perihelion. 
Bottom: The last 60,000 yr in more detail; substantial disintegra¬ 
tion is likely to take place in this orbital state. 

About half temporarily entered short-period orbits, and nine 
temporarily entered Earth-crossing orbits, often repeatedly. 
There were 65 Earth-crossing epochs in all, ranging from 
~100 to ~ 14000 yr, with mean duration ~2000 yr (Fig. [21). 
The dynamical evolution of one such clone, which included a 
short-period phase {P ~3.5-5.3 yr) of 50,000 year duration, 
is illustrated in Fig. [3] In these runs the bodies are treated 
as point masses and their physical evolution is neglected. If 
100 simulated Chirons, with mean lifetime 0.35 Myr, yield 65 
Earth-crossing episodes, then the frequency of such episodes 
is 65/100/0.35 Chiron“^Myr“^, giving recurrence time ~1.9 
Myr per Chiron clone. 

2.2 Estimates of the Centaur population 

Di Sisto & Brunini (2007) estimated the number of Centaurs 
larger than Chiron to be 360< N <650, of which there are 
4< N <7 within 18 au of the Sun (Centaurs are here defined 
as objects with perihelia between Jupiter and Neptune, ir¬ 
respective of semi-major axis). This is consistent with the 
discovery of four Centaurs larger than Chiron, three of them 
within 18 au of the Sun. If there are typically four objects 
of this size at any given time with similar orbital properties, 
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each with arrival probability into an Earth-crossing orbit of 
p ~0.5 Myr“^, then there is an expectation of one Chiron¬ 
sized Earth-crosser every half a million years or so. 

The size distribution of Centaurs may depend on the rel¬ 
ative contributions from their source populations. For long- 
period comets, Fernandez & Sosa (2012) find a cumulative 
radius R distribution N{R) oc R~‘‘, with s = 2.15 ± 0.75, 
consistent with s = 2 for an evolved population in the diam¬ 
eter range 20-200 km, expected from current models of plan- 
etesimal growth (e.g. Kenyon & Bromley 2012). Bauer et al 
(2013) analysed thermal infrared observations of 52 centaurs 
obtained from the Wide-held Infrared Survey Explorer and 
obtained a biased power law size distribution with index 
s = 2.7 ± 0.3. If we adopt a diameter 11=240 km (Stans¬ 
berry et al 2007), then with these distributions the number 
of comets in Chiron-like orbits, >100 km diameter, entering 
Earth-crossing orbits in a million years is ~6-14 times the 
Chiron rate of entry, or about one every 36,000-86,000 years. 
This hgure is likely to be conservative as additional con¬ 
tributions are expected from more remote Centaurs. Asher 
& Steel (1993), for example, followed the evolution of the 
large Centaur 5145 Pholus over 1 Myr using 27 clones, which 
yielded two particles entering Earth-crossing orbits for sev¬ 
eral tens of thousands of years, implying that a body in such 
an orbit has a 5-10% chance becoming Earth crossing within 
a million years. 

Dehning the system of Jupiter family comets (JFCs) 
as having orbital periods P <20 yr and Tisserand criterion 
T <3, Di Sisto et al (2009) estimate a population of A^=450 
with perihelia q <2.5 au. For a mass distribution n(m) dm oc 
dm in the radius range 1< R <10 km, the mass of the 
JFC system is then 

— 1, TUlo .... 

Mjfc = N{- -)- mi (1) 

2 — 0 ! mi 

where (mio,mi) are respectively the masses corresponding 
to JFCs of radius (10 km, 1 km) respectively. The size dis¬ 
tribution adopted by Di Sisto et al translates to a pop¬ 
ulation index o = 1.9. With this, and mio/mi = 1000, 
we find Mjfc ~8100mi. A comet of radius 1 km and 
density p = 0.5 gmcm“® has mass mi ~2.6xl0^^ gm. 
Hence Mjfc ~2.2x 10'^® gm. A similar exercise for the 
current near-Earth asteroids (NEA) yields, if they are 
largely composed of rock with density 2.5 gmcm“®, a mass 
~5xl0^® gm. A 250 km comet with density 0.5 gmcm“® 
has mass ^5x10^*^ gm, two thousand times that of the cur¬ 
rent JECs and a thousand times that of the NEAs. Thus in 
terms of Earth-crossing mass the present-day impact hazard 
represents, not a steady state, but a snapshot of a strongly 
flickering system. The fluctuation timescale is ~0.5 Myr for 
bodies of at least Chiron size, and ~0.03-0.1 Myr for comets 
at least 100 km in diameter. 


3 BREAKUP MODES 

The outgassing of an active Chiron-sized comet in say a 5 yr 
orbital period (Fig. [3]) is equivalent to that of 15,000 Halley¬ 
sized comets in Encke-like orbits, or in excess of 40 strong 
meteor showers per night, and the same by day. Apart from 
periods of dormancy, this will persist over the dynamical or 
physical residence time of the comet, whichever is shorter. 
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Figure 4. Part of Fig. [3] at higher resolution still, showing a 
transition from water snowline (sublimation at (;<2.5au) to a 
short-period orbit with g <1 au in less than 1000 years. 

Such intense atmospheric bombardment is expected to have 
consequences for the biosphere (Section|6]). This decay mode 
may be checked, albeit temporarily, by the development of 
a refractory porous crust. In the rubble mantle model of 
Jewitt (2002), the timescale for mantle growth for a water 
ice and carbon monoxide composition is of order a year at 
2.5 an and ~100 yr at 5 an. The transfer of a Centaur from 
an orbit with q = 2.5 au (the water snowline) to an Earth¬ 
crossing orbit may take place in less than 1000 years (Fig. [4]). 

‘Spontaneous’ splitting or fragmentation of very large 
comets are probably unimportant disintegration routes since 
their escape velocities (~70 ms“^ for a 250 km comet, ~30 
ms”*^ for a 100 km one) are in excess of the likely separation 
speeds of fragments. Disintegration is more likely to take 
place through sublimation, or tidal disruption by the Sun or 
Jupiter. The simplest breakup mode, namely generation of 
an intense meteoroid stream along the track of the comet, 
has been discussed by Clube & Napier (1984). According to 
Opik (1951), a comet loses 4/y^ m of surface per revolution, 
a figure not very different from that found empirically by 
Di Sisto et al (2009). This yields an output of 10^® gm of 
gas and dust per revolution for a Chiron-sized body in an 
Encke-like orbit, but could rapidly produce a dormant stage, 
so preserving the integrity of the comet as its perihelion is 
driven downwards (e.g. Rickman et al 1990). 

Tidal disruption by Jupiter is also a common pro¬ 
cess, several examples having occurred in recent times (e.g. 
D/1993 F2 Shoemaker-Levy 9). The timescale for collision 
of a 1 km JFC on Jupiter is about a decade (Zahnle et 
al. 1998, Sanchez-Lavega 2010), implying that passage of a 
comet within its Roche limit (~3.4i?j) is virtually an annual 
event. It is thus likely that the JFC system is replenished 
sporadically by the disintegration of large comets passing 
within the Jovian Roche lobe (Bredichin 1889). Pittich & 
Rickman (1994) have shown that fragments of a JFC so dis¬ 
rupted would lose their memory of a common origin within 
a few thousand years. A ~10® enhancement of the JFC pop¬ 
ulation by the disruption of a 250 km comet would feed into 
a corresponding increase in the NEO cometary population 
over the physico-dynamical decay time of the JFC. 

A common fate of comets thrown into short-period or¬ 


bits, if not disintegration through sublimation, is tidal dis¬ 
ruption by the Sun (e.g. the ~100 km Kreutz progenitor). 
Destruction through tides is a volume rather than a surface 
phenomenon, and is only weakly dependent on the size of 
the comet. Bailey et al (1992) have shown that, under the 
influence of slowly acting secular perturbations of the giant 
planets, ~5-15% of comets with initial perihelion distances 
q <2 au may evolve into a sungrazing or sun-colliding state. 
This evolution occurs at intervals of order 10® revolutions, 
about a tenth of the dynamical ejection timescale. The tidal 
disruption of a comet entering the Roche limit, ~3.57?0, may 
create a hierarchy of fragmentations. Passage through such 
clusters, before they disperse into the background, will gen¬ 
erate one or more series of sharp bombardments, analogous 
to meteor storms but of much greater intensity, throughout 
the lifetime of the comet. 


3.1 Encke-like orbits 

Ipatov & Mather (2004), following the orbital evolution of 
26,000 fictitious Jupiter-crossing objects, found that bodies 
in orbits close to that of Comet 2P/Encke had mean collision 
probabilities with the terrestrial planets about two orders of 
magnitude greater than those of other comets, remaining in 
typical NEO orbits for millions of years. Levison et al (2006), 
following orbits from the Kuiper belt, deduced that at any 
one time there should be one to a dozen objects in Encke-like 
orbits, neglecting both physical decay and nongravitational 
effects. They found that the time taken for a JFC to attain 
an Encke-like orbit, dynamically uncoupled from Jupiter, is 
0.1 Myr, two orders of magnitude longer than its physical 
lifetime, but proposed that 2P/Encke might have survived 
destruction through a period of dormancy. 

Unless dormancy is 100%, however, non-gravitational 
forces may operate and enhance the rate of transfer sub¬ 
stantially in combination with mean motion resonances with 
Jupiter (Steel & Asher 1996; Fernandez et al 2004; Pittich 
et al 2004). The acceleration due to such forces varies as 
area/mass, i.e. /B? oc 1/i?, but the fuel available oc R® 
and unless permanent dormancy intervenes, long term non¬ 
gravitational perturbations (acceleration times time) are 
thus more effective in perturbing large comets (oc R^). In the 
event of a perihelion disintegration, for example, the cross- 
sectional area of fragments reaching the Earth’s distance is 
typically about 10,000 times that of the comet itself, for an 
initial dispersion speed of 5 ms“^ (Napier 2010). However 
the transition probability from the JF to Encke domains 
remains uncertain to order of magnitude. 

In summary, on the arrival of a giant comet in an Earth¬ 
crossing orbit, whether Jupiter family or Encke-like, sev¬ 
eral modes of disintegration are possible which may act to 
enhance the terrestrial bombardment rate by two or three 
powers of ten over background, during its physico-dynamical 
lifetime. Direct impact of an incoming large comet is much 
less probable than interaction with its debris, which may 
include many short-lived fragments, and a correspondingly 
spiky distribution of terrestrial impactors is expected. This 
can be tested by examining the impact cratering record. 
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D (km) 

age (Myr) 

crater 

no. 

epoch (Myr) 

30-120 

0.78±0.1 

Austr. tektites 

1 


14 

0.9±0.1 

Zhamanshin 

1 


10.5 

1.07±0.0 

Bosumtwi 

1 


3.4 

1.4±0.1 

New Quebec 

1 

l.li0.03 

52 

2.5±2.5 

Kara-Kul 



18 

3.5±0.5 

El’gygytgyn 



10 

5±1 

Karla 



8 

5±3 

Bigach 



24 

15.1±0.1 

Ries 



40 

35.3±0.1 

Chesapeake 

2 


90 

35.7±0.2 

Popigai 

2 


28 

36.4±4 

Mistastin 

2 


7.5 

37.2±1.2 

Wanapitei 

2 

35.4i0.2 

15 

42.3±1.1 

Logoisk 



9 

46±3 

Ragozinka 

3 


25 

49.0±0.2 

Kamensk 

3 


45 

50.50±0.76 

Montagnais 

3 

49.1i0.2 

12.7 

58±2 

Marquez 



150 

64.98±0.05 

Chicxulub 

4 


24 

65.17±0.64 

Boltysh 

4 

65.0i0.05 

65 

70.3±2.2 

Kara 



35 

74.1±0.1 

Manson 



6.5 

81.Oil.5 

Wetumpka 



19 

89.0i2.7 

Dellen 



13 

99i4 

Deep Bay 



9 

121.0i2.3 

Mien 

5 


12 

123il.4 

Vargeao Dome 

5 

122.5il.2 

55 

128i5 

Tookoonoka 



40 

142.0i2.6 

Mjolnir 

6 


22 

142.5i0.8 

Gosses Bluff 

6 


70 

145.0i0.8 

Morokweng 

6 

143.7i0.6 

3.2 

165i5 

Zapadnay 

7 


40 

167i3 

Puchezh-Kat. 

7 

166.5i2.6 

40-50 

200 il 

seismite 

8 


23 

201 i 2 

Rochechouart 

8 

200.2i0.9 

85 

214il 

Manicouagan 



40 

254.7il.3 

Araguainha 



52 

376.8il.7 

Siljan 

9 


8.5 

378i5 

Ilyinets 

9 


15 

380i5 

Kaluga 

9 

377.2il.5 

6 

445i2 

Pilot 




Table 1. Impact structures, selected as discussed in the text. 
Data from the Earth Impact database, with the addition of the 
Australasian tektites and a seismite impact marker at 200 Myr. 
Possible impact episodes (no.) are listed, epochs being given as 
the weighted mean age of the associated craters, with an unbiased 
weighted estimate of the sample standard error. Updated from 
Napier (2006); cf Napier (2014). 


4 BOMBARDMENT EPISODES 

Of the 184 confirmed craters in the Earth Impact Database, 
accessed in 2014 September, 39 were culled with D >3 km, 
dating accuracy 2cr <5 Myr, age <500 Myr. Two non-crater 
impact structures were added: (a) Simms (2003, 2007) has 
identified a seismite layer which appears to have been cre¬ 
ated from another large impact, yielding perhaps a 40-50 
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Figure 5. Age distribution of Table [T] impact structures. Data 
have been smoothed by a rectangular window 6 Myr wide, in 
1 Myr steps. Estimates of the relative rate of decline as one goes 
into the past are shown by the curves, derived as described in the 
text; the upper curve is the Sturges fit, the lower is the Rice one. 

km unidentified crater, at least 600 km west or northwest of 
central Britain; and (b) the Australasian tektites have been 
added; their source crater too has not been identified, but 
published estimates of its size range from 30 to 120 km (Lee 
& Wei 2000). Thus the database comprises 41 well-dated 
impact structures covering the last 500 Myr. This is listed 
in Table[T]and shown in Fig. [5] The solid lines are estimates 
of the decline of observed crater numbers with time, and 
are cubic fits to histograms of the database with bin widths 
respectively k = 7Q Myr and 100 Myr. There is no unique op¬ 
timum histogram giving a broad trend; the former is optimal 
under the ‘Rice rule’ k = , n the number of data (Lane 

2014), and the latter derives from the rule of Sturges (1926), 
namely k = log 2 n -|- 1. A somewhat larger dataset (n = 50 
craters with ct <10 Myr)) was taken in obtaining these sec¬ 
ular trends. Discovery completeness declines to 50% within 
120 Myr for impact craters with diameters D >40 km, and 
within 70 Myr for smaller craters. Identification of impact 
craters with D >20 km is only ~15% complete over the 500 
million year period, declining to ~5% complete for impact 
craters >10 km across. Nearly all such craters have been 
discovered on land. Evidently the record of impact craters 
is erased by erosion, sedimentation and subduction. In Ta¬ 
ble [T] we define the epoch of an episode by the weighted 
average date of the impact craters belonging to it. The in¬ 
trinsic duration of each bombardment episode could be as 
little as a few millennia or as much as 1 or 2 Myr if caused by 
the breakup of a single, large Jupiter family or Encke-type 
comet, for example via the Jovian Roche lobe. 

Fig. [6] is a plot of the probability density distribution 
of the listed impact structures, using Gaussian kernels cor¬ 
responding to the individual age errors. The Figure reveals 
that there is a tendency for impact craters to cluster in time. 
Fig. [7] shows the differential ages of the real dataset along 
with the average of 1000 equivalent datasets obtained by 
random extraction from the erosion curves (which of the two 
curves of Fig. [S] is employed makes no appreciable difference 
to the outcome). A slight excess of crater pairs relative to a 
random distribution is evident, and persists throughout the 
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age (Myr) 

Figure 6. Probability density distribution of impact structures 
listed in TablJTI using Gaussian kernel estimators corresponding 
to individual standard errors. Top: 0-100 Ma; bottom: 100-500 
Ma. The putative multiple impact episodes of TablcJJare labelled, 
with the corresponding number of impact structures in brackets. 


timespan of the dataset. Table [T] reveals that 26 of the 41 
impact structures listed belong to about nine temporal clus¬ 
ters with 5t<2 Myr, each of whose members are generally of 
identical age to within their standard deviations. Rejecting 
impact craters ages over 4 Myr from the mid-point of the 
nearest cluster yielded the weighted mean cluster ages listed 
in Table 1. The maximum deviation from the weighted mean 
within each episode averages 1.6 Myr. 

There are 27 pairs of craters with age separations 6t < 
2 Myr in the dataset (a triplet of craters bunched within 
2 Myr would yield three pairs). To test whether this pair 
count is compatible with the random arrival of single im- 
pactors, comparison was made with expectations from syn¬ 
thetic data. Two models were based on 10,000 extractions 
in batches of 41 from each of the secularly declining curves 
of Fig. [S] the number of crater pairs being counted in each 
dataset so constructed. A further approach was to apply 
random perturbations to the crater ages, enough to wipe 
out any clustering but not enough to significantly change 
the overall distribution. Ages were varied randomly by up 
to ±50 Myr subject to keeping the data within (0, 500) Myr. 
The results are shown in Fig.0. The number of trials where 



dt (Myr) 



dt (Myr) 



dt (Myr) 


Figure 7. Distribution of differential ages for the 41 impact 
craters (histogram), and for 1000 synthetic datasets extracted 
from the smoothed (Rice) distribution, normalised to the real 
dataset. Top to bottom: craters in the age ranges 0-35, 35-500 and 
100-500 Myr respectively. An excess of clustering up to 2 Myr (rel¬ 
ative to a random distribution) is found throughout the record. 
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number of pairs 



Figure 8 . Three simulations, each involving 10000 trials, in which 
pairs of craters are counted having age difference 5t <2 Myr. 
From highest to lowest, the curves correspond to models in which 

(a) the actual crater ages are randomised by up to ±50 Myr; 

(b) erosion is fitted to bins 100 Myr wide; and (c) fitted to bins 
70 Myr wide. The actual dataset (arrow) has 27 such pairs. The 
synthetic datasets all yield similar probabilities ~ 2 X 10“® that 
this number would be obtained by chance. 


Figure 9. A series of runs, in each of which 100,000 synthetic 
datasets are created from the real one by applying random per¬ 
turbations up to a maximum (‘fuzz’) to individual crater ages. 
For each 100,000-trial run, the number n of trials yielding at least 
27 pairs (<5t <2 Myr) is recorded. As the time resolution of the 
record is degraded, it asymptotically approaches a random con¬ 
figuration which reproduces the observed number of crater pairs 
in only 0 . 2 % of trials. 


>27 pairs out of 10,000 were obtained were, for the Rice 
and Sturges secular fits, and for the randomized datasets, 
respectively (37, 57, 42). Thus broad secular trends make 
little difference to the counting of closely separated crater 
pairs, all three yielding departures from randomness with 
probability ~99.8%. 

In creating synthetic datasets by applying random per¬ 
turbations </ Myr to the real one, any significant difference 
which emerges can only be due to structures on timescales 
</, since by construction the datasets are otherwise identi¬ 
cal. In the limit / —0, synthetic and actual datasets become 
identical and no signal can be detected. By progressively de¬ 
grading the resolution (increasing /), memory of the struc¬ 
ture is progressively lost from the synthetic datasets, and 
an upper limit can be placed on the timescale of whatever 
non-random structure is present. Fig. [5] shows the outcome 
of this exercise for sets of 100,000 trials in each of which the 
crater ages were randomly varied by up to / Myr and the 
number of pairs {St <2 Myr) was counted. Randomisations 
of ±20 Myr or more in the crater ages yield the observed 
number of pairs with probability 0.2%. However, it can be 
seen that even on timescales < 10 Myr, nonrandom struc¬ 
ture is present with probability ~99%. For the 16 craters 
older than 100 Myr listed in Table [T] the mean separation is 
~25 Myr and yet there are 6 crater pairs with St <1.6 Myr. 
Randomization trials with these data yielded 6 or more such 
pair counts in ~0.5% of 10,000 cases. 

Comparison was then made with models comprised of 
bombardment episodes, generated under the following as¬ 
sumptions: 

(a) Large comets were assumed to arrive into Earth-crossing 
orbits (JFC or Encke) at random, and disintegrate. 

(b) Impact craters were randomly generated from the frag¬ 
ments in numbers proportional to the comet mass, which 
was extracted randomly from a population index a ~ 1.9 
appropriate to the JFC population, to some maximum. This 
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number of pairs 

Figure 10. Four simulations, each involving 10000 trials, in which 
pairs of craters are counted having age difference St <2 Myr. The 
actual dataset has 27 such pairs. The curve on the left shows 
the distribution of such pairs when clustering is purely random. 
The others correspond to bombardment episodes of maximum 
amplitudes 10, 20 and 50, structured as described in the text. 


latter arises from a combination of physical and dynamical 
evolution (Di Sisto et al 2009). For models in which im- 
pactors are the disintegration products of cascading frag¬ 
mentation, a is probably somewhat smaller (Sekanina & 
Chodas 2007; Ishiguro et al. 2009). 

(c) The measured duration of a bombardment episode was 
taken to be 2 Myr, determined largely by the dating uncer¬ 
tainties (Table 0. 

(d) Typically, in each trial, 200 craters were generated in 
this way over a 500 Myr period. They were then reduced 
in numbers randomly, in proportions matching the observed 
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decline of detected craters derived from the Rice cubic fit, 
until 41 remained. 

Fig. [To] shows the outcome of trials in each of which 
10,000 fictitious datasets were constructed. Each trial had, 
successively, a maximum possible amplitude ric = 2, 3, 4 
... for episodes of bombardment. It is clear from the Figure 
that the observed count of 27 pairs with 5t < 2 Myr is rare 
under the hypothesis Ho of no clustering, but is comfortably 
encompassed by a wide range of bombardment amplitudes. 
The probability of no impact clustering, given the dataset 
D, can be expressed as 


P{Ho\D) 


P(D\Ho)Pprior{Ho) 

PiD\Ho)Pprior{Ho) + PiD\Hl)Pprior{Hl) ^ ^ 


where 


P{D\Hl)PprioriHl) = j:P{D\H')Pprior{H') (3) 

Hi represents the clustering hypothesis, Pprior{Ho) and 
Pprior{H') represent prior probabilities in Bayesian termi¬ 
nology and the summation is taken over all viable cluster 
models H' corresponding to = 2, 3, 4 .. .. 

We adopt the prior hypothesis that the data are equally 
likely to be clustered or non-clustered. The relative probabil¬ 
ities are calculated from the occupancy of the n = 27 bin un¬ 
der the various hypotheses. On the no-clustering hypothesis, 
the bin with 27 hits has occupancy 15±5 ‘hits’ out of 10,000 
trials: P{D\Ho) = 1.5 x 10“®. The summation of the bin 27 
populations for ric = 2, 3, ... 100 yields a total occupancy of 
~15,000 for a million trials and hence P{D\Hi) = 1.5 x 10“^. 
Thus the posterior odds ratio P{Hi\D)/P{Ho\D)'^10, a pos¬ 
itive result favouring the existence of clustered impacts. 


4.1 Bombardment intensities and causes 

Although the existence of impact episodes is thereby indi¬ 
cated, their intensity is not well constrained. The widths of 
the bombardment episodes in the trials, typically 2 Myr, 
may reflect dating uncertainties rather than their true du¬ 
ration. An amplitude of factor three above background over 
2 Myr translates to a factor of 300 over background for a 
comet disintegration of 20,000 years’ duration, or 3000 over 
2000 years. The only impact episode for which we have pre¬ 
cise relative dating is that of 65 Myr BP. The Chicxulub 
structure of that era is the largest impact crater known to 
have been formed over the Phanerozoic, and is generally as¬ 
sociated with the end-Cretaceous extinctions with estimated 
species loss of ~57-83%. Although originally ascribed to an 
asteroid impact, iridium and energy constraints yield a best 
estimate that the impactor was a comet 10-80 km in diame¬ 
ter (Durand-Manterola & Codero-Tercero 2014). However a 
detailed stratigraphic study of the Ukrainian Boltysh crater 
by Jolley et al (2010) has shown that it was formed only 
2000-2500 years before the Chicxulub impact in the Gulf of 
Mexico (the radiometric age given in Tabl^H has standard 
deviation ±0.6 Myr). This short timescale is consistent with 
the breakup of a large short-period comet yielding N Earth¬ 
crossing fragments with dynamical or physical lifetimes L yr 
such that NL>C, where C represents the mean lifetime of 
a fragment before collision with the Earth. For a comet in 
an Encke-like orbit, (7^200 Myr and 10® km-sized fragments 
lasting 2000 yr (with total mass corresponding to a nominal 


diam (km) 

in 

out 

> 20 

15 

7 

10-20 

4 

6 

3-10 

6 

3 


Table 2. Membership of bombardment episodes listed in Table]!] 
There is no evidence of a systematic trend with crater diameter. 


70 km diameter comet) or 10^ fragments lasting 20,000 yr 
(a 30 km comet) would suffice. 

Other contenders for bombardment episodes are comet 
showers from the Oort cloud, or impacts from the disinte¬ 
gration of asteroids in the main belt, but neither mechanism 
seems capable of reproducing the sharpness and amplitude 
of the observed episodes. Kaib & Quinn (2009) constrain 
the strength of comet showers through stellar penetrations 
of the inner Oort cloud, by noting that too high a comet 
number density in the inner Oort cloud would overproduce 
long period comets. They find that there has been no more 
than about one shower of Eocene-Oligocene strength over 
the past 500 Myr. Terrestrial impacts from the breakup of 
main belt asteroids are too infrequent by order of magnitude 
to account for the episodes (Napier 2006). The main aster¬ 
oid families generally yield 0-4 terrestrial impactors >1 km 
across, arriving over 5-90 Myr (Zappala 1998). 

The question arises whether the tendency to clustering 
is a function of crater diameter. In Table [5] the craters are 
partitioned by size and by membership or otherwise of an 
impact cluster. No significant correlation is evident between 
crater size and cluster membership: the tendency for im¬ 
pactors to arrive in groups extends down to sub-kilometre 
bodies. 


4.2 The most recent episode 

The apparent spike in the bombardment rate of the past few 
million years could be seen as due simply to the relative com¬ 
pleteness of discovery of recent craters. However, between 
the eight well-dated craters < 5 Myr in age and the Chesa- 
peake/Popigai group at ~36 Myr, there is only one large im¬ 
pact structure, the 24 km Nordlinger Ries at 15 Myr (from 
an achondritic impactor). The Earth Impact Database, with 
184 craters, yields only one additional crater, the 3.6 km di¬ 
ameter Colonia with an age in the range 5< t <35 Myr. On 
a random impact scenario, the lunar cratering record indi¬ 
cates that bodies large enough to make 20 km craters strike 
the Earth every 600,000 years (Bland 2005). Allowing for 
ocean impacts, there is then an expectation that ~17 land 
craters of at least this size should have been formed over the 
30 Myr interval (Bland 2005), most of which should have 
been discovered since this a geologically young period. The 
gap therefore appears to be real. 


5 COMPARISON WITH MARINE 
EXTINCTION RATES 

For comparison with extinction data, the 500 Myr marine 
fossil record database assembled by Bambach (2006) was em¬ 
ployed. This came from the Paleobiology Database and the 
compendium by Jack Sepkoski: discussion and references are 
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Figure 11. Bombardment episodes (asterisks) and relative ex¬ 
tinction intensities of marine genera from Bambach (2006). Up¬ 
dated from Napier (2014). 


No. extinction 

episode 

At 

boundary 

1 

1.8 

1.1 

0.7 


2 

33.9 

35.4 

1.5 

end Eocene 

3 

65.5 

65.0 

0.5 

end Cretaceous 

4 

93.5 

- 

4.2 


5 

145.5 

143.7 

1.8 

end Jurassic 

6 

183.0 

- 

- 


7 

199.6 

200.2 

0.6 

end Triassic 

8 

251.0 

254.7 

3.7 

end Permian 

9 

260.4 

- 

- 


10 

322.2 

- 

- 


11 

359.2 

- 

- 


12 

374.5 

377.2 

2.7 

end Devonian 

13 

385.3 

- 

- 


14 

443.7 

445.2 

1.5 


15 

445.5 

445.2 

0.3 

end Ordovician 

16 

488.3 

- 

- 


17 

494.7 

- 

- 


18 

498.9 

- 

- 



Table 3. Epochs (Myr ago) of marine genera extinction events 
considered by Melott & Bambach (2010) to be genuine mass ex¬ 
tinctions, along with those impact episodes of Fig. [ 6 ] occurring 
within At <5 Myr of the extinctions. Seven of the twelve such 
events occur at major geological boundaries. 

found in Melott & Bambach (2010). The impact crater ages 
employed have been determined radiometrically, whereas the 
marine extinctions are generally based on stratigraphy and 
contain rounding errors due to the assignment of extinctions 
to the nearest geological substage. The relative extinction 
intensities of genera are plotted in Fig. 1111 along with the 
bombardment episodes. According to Bambach, 18 events in 
the last 500 Myr may properly be described as mass extinc¬ 
tions. In Table[3]these extinction peaks are listed along with 
those bombardments which occur within ±5 Myr of a peak. 
The well-dated single impacts in Fig. [6]are included. Two of 
the mass extinctions, events 14 and 15, are extremely close 
together and share a bombardment episode. For statistical 
purposes they are here counted as one. If we regard each 
bombardment episode as a single event and include the sin- 
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gle impacts, there are 14 impact episodes, seven of which lie 
within ±3 Myr of an extinction episode, six within ±2 Myr 
and four within ±1 Myr. Seven coincide with the termina¬ 
tions of major geological boundaries to within the dating 
errors, namely the Ordovician, Devonian, Permian, Trias- 
sic, Jurassic, Cretaceous and Eocene. The bombardments 
comprise 11 multiple episodes (Table [1]) and three singles. 

Monte Carlo trials were carried out in which 14 syn¬ 
thetic bombardment episodes were randomly extracted over 
500 million years, and the number which lay within At Myr 
of the Bambach mass extinction (Table[3]) was recorded. Not 
more than one episode could be assigned to an extinction, 
since they would not be distinguishable. In reality, uncer¬ 
tainty attaches also to the precise dating of the extinctions, 
and this may be soaked up in At. Trials for each of the 
thresholds At = ±(3,2,1) Myr yielded probabilities (0.032, 
0.005, 0.012) that the coincidences are due to chance. Ex¬ 
cluding the first 30 Myr of the extinction and bombard¬ 
ment records (in case the recent ‘episode’ is actually discov¬ 
ery bias), makes no essential difference to the probabilities. 
These coincidences independently validate the reality of the 
cratering episodes since, if they were statistical artefacts, 
there would be no reason to expect them to correlate with 
the observed mass extinctions. They also demonstrate that 
bombardment episodes have been major factors in the evo¬ 
lution of life since the pre-Cambrian explosion. 

There is continuing debate amongst palaeontologists 
about whether the Cretaceous-Palaeogene extinctions were 
sudden and caused by the Chicxulub impact, or whether 
there was a more continuous decline over < 1 Myr of species 
caused by global cooling, sealevel fluctuations and the main 
phase of the Deccan vulcanisms (MacLeod 2013). No other 
mass extinction has yet been identified with a single large 
crater, and most or all of the impact peaks are due to im- 
pactors too small to have caused major mass extinctions 
on their own. For example, the Araguainha crater in cen¬ 
tral Brazil is only 40 km across and was itself incapable 
of causing the near-extinction of life at the Permo-Triassic 
boundary 250 million years ago; the end-Ordovician peak at 
445 Myr represents two craters only a few kilometres in di¬ 
ameter; the Triassic peak at 200 Myr represents two impact 
structures with diameters in the probable range 20-40 km, 
and so on. The collectively strong association of such minor 
impactors with major extinctions implies that they are prox¬ 
ies rather than individual causes of the mass extinctions. 


6 TERRESTRIAL CONSEQUENCES 

The stepwise nature of several mass extinctions has often 
been noted, and the hypothesis of multiple impacts as their 
cause, based on purely biostratigraphic considerations, is 
not new (e.g. the Late Devonian crisis of 365-368 Myr BP: 
McGhee 1996). Some of these ancient events were highly de¬ 
structive, with intensities that came close to extinguishing 
life on Earth. There is, however, a spectrum of intensities, 
as there is a spectrum of comet masses, and the correla¬ 
tion evident in Fig. [TT] and Table |3] suggests that episodes 
of bombardment may exert a more continuous control over 
the evolution of species than would be expected simply from 
rare, isolated impacts. A progressive ‘toughening’ of marine 
genera with time is evident from the slope of the curve in 
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Fig. 111! There is no obvious astrophysical explanation for 
this decline in extinction rate, but no need for one as several 
palaeontological hypotheses have been advanced to account 
for it (MacLeod 2013). 

Birks et al (2007) have pointed out that the most dam¬ 
age which can be done to the biosphere with the least impact 
energy is to deplete the stratospheric ozone layer, so increas¬ 
ing the level of UVB radiation reaching the surface of the 
Earth. Depletion by way of impact occurs partly through the 
activation of halogens thrown into the stratosphere through 
salt-water impact, and partly through the injection of NO 
created in the impact fireball. They consider that the impact 
of a ~0.5 km asteroid would have serious implications for the 
biosphere, and that such impacts occur at ~0.1 Myr inter¬ 
vals. The effects of enhanced UVB on terrestrial organisms 
have been discussed by numerous authors. Melott & Thomas 
(2011) have pointed out that a current anthropogenic 3% de¬ 
pletion in global stratospheric ozone is causing measurable 
damage to the biosphere (Hader et al 2007). They consider 
that a 30% depletion is well above the lethality level for 
phytoplankton, which lie at the base of the oceanic food 
chain. Such a depletion would cause a major marine extinc¬ 
tion event. Melott et al (2010) have also pointed out that 
cometary airbursts may create nitrous oxides in the strato¬ 
sphere, depleting ozone. In the present scenario, bombard¬ 
ment episodes may enhance the frequency of small impacts 
by two or three powers of 10 over background during the 
disintegration of a giant short-period comet. 

A further possible effect is climatic perturbation due to 
the dust production expected in the course of disintegration. 
If cometary fireballs are reduced to micron-sized particles 
on entering the atmosphere (Klekociuk et al 2005), then a 
current infall of cosmic dust of 40,000±20,000 tons a year 
(Love & Brownlee 1993) yields an atmospheric optical depth 
r~8xl0~^ for an absorption coefficient «;~10“® cm^ g“^. 
This figure may rise to r~0.2 over the course of a giant 
comet disintegration (Napier 2001), raising the prospect of 
repeated pulses of global cooling each of millennia in dura¬ 
tion. 


7 SUMMARY AND CONCLUSIONS 

As with the periodicity issue, the analysis of bombardment 
episodes is inhibited by small number statistics. Neverthe¬ 
less positive evidence for such episodes appears in the impact 
cratering record. They are found to be tightly correlated 
with substantial marine extinctions at the level of genera. 
To date only the Cretaceous-Tertiary extinctions have been 
widely identified as being caused by a single large impact, 
and even that claim is disputed by many palaeontologists 
(Keller et al 2012, MacLeod 2013, Keller & Kerr 2014). Sev¬ 
eral large impact craters, which might have been expected 
to have caused mass extinctions, are associated with only 
minor, regional ones. The main asteroid belt seems inad¬ 
equate to supply the larger impactors, and large cometary 
populations have now been revealed extending to beyond the 
planetary system, with the potential to be thrown into the 
near-Earth environment. The top-heavy mass distribution of 
these comets, their likely disintegration history when thrown 
into short-period orbits, and the potential of their debris to 
cause deleterious effects, make them likely candidates for 


the mass extinctions observed in the terrestrial record and 
the associated geological disturbances. The mechanisms in¬ 
volved are a matter for detailed stratigraphic studies; it is 
suggested here that prolonged atmospheric perturbations 
arising from fireball storms and dusting are the most energy- 
efficient means of collapsing food chains, yielding both ma¬ 
rine and land extinctions. 
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